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Abstract I discuss the properties of electron states in 
amorphous Si based on large scale calculations with real- 
istic several thousand atom models. A relatively simple 
model for the localized to extended (Anderson) transition 
is reviewed. Then, the effect of thermal disorder on local- 
ized electron states is considered. It is found that under 
readily accessible conditions, localized (midgap or band 
tail) states and their conjugate energies may fluctuate 
dramatically. The possible importance of non-adiabatic 
atomic dynamics to doped or photo-excited systems is 
briefly discussed. 

1. Introduction 

There have been a very large number of studies on 
electron states in amorphous materials, and especially 
amorphous silicon [1]. These computations have ranged 
from simplified tight-binding approaches to first princi- 
ples methods using density functional theory in the lo- 
cal spin density approximation (LSDA) [2]. A difficulty 
with many calculations of electron states in amorphous 
Si is the usual challenge of obtaining a realistic starting 
model. Amorphous Si is well known to be almost entirely 
four-coordinated with fewer than one atom in 10'^ be- 
ing non-ideally coordinated. However, few if any models 
fabricated entirely with molecular dynamics (MD) (in- 
cluding first principles MD) have a defect concentration 
at all consistent with experiment. In particular, models 
obtained by a quench from the melt are particularly un- 
realistic, often with 10% or more coordination defects in 
a cell of order 100 atoms. Such a model does not posess 
a true optical gap, but rather just a deep minimum in 
the density of states near the Fermi level. This result is 
inconsistent with electron spin resonance (ESR) measure- 
ments and most other estimates of defect concentration in 
the material. Such models are probably best interpreted 
as poor metals rather than semiconductors, since there 
is so much banding between impurities that the defect 
states are fully extended through the supercell. The ori- 
gin of this trouble is easy to ascertain: the liquid (a pre- 
dominantly sixfold metal) is topologically quite distinct 
from amorphous (or cystalline) Si, and any rapid quench 
is certain to freeze in too many remnants of the liquid 
(such as overcoordination defects, which are more com- 
mon than dangling bonds in most published MD simula- 
tions [2]). Because of the qualitative difference in liquid 
and amorphous topology, it is perhaps not surprising that 
a-Si cannot be formed experimentally from a quench from 
the liquid. We have shown elsewhere [2] that a proper de- 
scription of defect wave functions (with localization con- 



sistent with estimates from ESR) requires both a realistic 
model (with a very small defect concentration less that 
1 % ) and LSDA, to estimate the spin localization of an 
isolated defect. Fortunately realistic models (meaning 
in uniform agreement with structural, vibrational and 
electronic/optical experiments) do exist, but are made 
using Wooten-Weaire- Winer [3] methods. An excellent 
set of these models (with 216, 512 and 4096 atoms) was 
produced by Djordjevic, Thorpe and Wooten [4]; we use 
these in the studies presented here. 

It is now established that the midgap defect state is 
due to danghng bonds. The floating bond conjecture 
(that five-fold coordinated atoms play an important role 
in a-Si) has not been disproven, though any such states 
are shifted well toward the conduction band edge [2] and 
are not nearly so localized as the sp^ dangling bonds. If 
such floating bond defects exist, they might form reso- 
nances with the conduction tail states, which could have 
interesting implications for transport in n-type materials. 

Less progress has been made on band tail states, be- 
cause a good description of these requires models with 
large numbers (several hundred to thousands) of atoms. 
To obtain a qualitative understanding of the Anderson 
transition in a three-dimensional network we used a 4096 
atom supercell model [4] (a cube of amorphous Si about 
4 nm on a side). The band tail states have structure on a 
length scale several times the nearest neighbor distance 
(depending on the energy), and an enormous model (by 
the standards of electronic structure calculations) is es- 
sential to disentangle the "simple" features of the band 
tail tates in amorphous Si. In section 3 we describe the 
properties of these states. 

Finally we discuss the electronic effects due to ther- 
mal disorder superposd on a topologically disordered net- 
work [5] . We observe from a straightforward thermal sim- 
ulation of an amorphous Si network that well-localized 
states can have their energies modulated by tenths of an 
eV (interestingly, an order of magnitude larger than the 
thermal energy itself). Also, the states themselves can 
fluctuate spatially as we describe in section 4 below. This 
feature is of importance to transport, optical and doping 
studies. 

2. Methodology 

For the calculations of electronic tail states in the 4096 
atom model, we used a tight-binding Hamiltonian and 
a sophisticated Lanzcos scheme to compute about 500 
states near the optical gap [6] (the dimension of the 
Hamiltonian matrix is 16384). A technical point of some 
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importance is that a "simple" Lanczos scheme designed 
for computing states near a selected energy, such as the 
folded spectrum method (computing states from a shifted 
and squared H matrix) is not adequate to compute more 
than about 20 states here because the spectrum becomes 
exceedingly dense for energies out of the middle of the 
gap. What is effective is a "shift and invert" strategy, 
which spreads the spectrum out effectively, enabling the 
reliable calciilation of many more states. Thus, one com- 
putes the matrix G{E) = {El — H)~^; and performs 
an elementary Lanczos calculation on the matrix G (to 
obtain spectral information near E; H is the Hamilto- 
nian matrix, 1 is the unit matrix). It is amusing that 
this method was apparently pioneered by engineers [7] 
instead of physicists, since the object G is so familiar 
to workers in condensed matter physics (it is the elec- 
tronic Greens matrix!). It is not a priori obvious that 
the expense of inverting {El — H) is worthwhile, but be- 
cause H is sparse this is not prohibitively expensive, and 
the transformed spectrum is well conditioned for Lanczos 
(the eigenvalues of G(£^) are well separated near E even 
if the spectrum of H is dense near E). 

For dynamical simulation, we use the methods of 
Sankey and coworkers [8] "Fireball96" , which is discussed 
in the context of amorphous Si elsewhere [5]. 

3. Model of the Anderson transition in amor- 
phous Si 

I have discussed the details of the computation of gap 
and band tail states in the 4096 atom model of amor- 
phous Si elsewhere [6]; here I dwell upon the "physics" 
of the calculation. In a nutshell, we found that (1) states 
near the middle of the gap were exponentially local- 
ized (as one would expect); (2) for energies shifted from 
midgap toward the valence tail, there was a tendency to 
find system eigenstates with rather similar charge den- 
sities, and consisting approximately of two separated 
"blobs" of charge; (3) deeper into the valence tail we ob- 
served that the states still were discernibly "built" from 
clumps or clusters of charge, with a clear separation be- 
tween these clumps; (4) deep into the valence tail we 
saw states that were well extended, essentially uniformly, 
through the cell. See Figs. 1 and 2 in Ref. [6]; (5) The 
states are very nonlocal away from midgap, in contradic- 
tion with the defect pool model; it is clear that this must 
happen since the states in the bands are extended; this 
computation shows that states are really well localized 
only within of order ^ 0.6 — 0.7eV from midgap. 

The preceding discussion suggests that an alternate 
basis or representation may be worthwhile to describe 
the electron states of a-Si, especially near the gap. As 
usual in quantum mechanics, one is at liberty to select 
any complete (or at least approximately complete) set 
to represent the quantum mechanical states. The choice 
is normally made on grounds of efficiency: basis func- 
tions are selected which captiirc much of the character of 
the eigenfunctions, so that only a small number of basis 



functions is required to faithfully represent the eigenfunc- 
tions. The classic example is the method of orthogonal- 
ized plane waves in band theory [9]. Likewise here, since 
the eigenvectors in the gap and near the tails are clearly 
composed of "lumps" or clusters of charge, a basis con- 
sisting of such cluster states is appcialing. 

I take the view that the cluster states are localized 
cigeufunctions which arise from idealized isolated defects. 
That is, one can imagine all possible distortions of a con- 
tinuous random network. Select one such distortion (say 
a strained bond angle removed from the distribution of 
bond angles in the system). Suppose further that the 
strain is large enough to make the state bound (eg. lo- 
calized) with exponential tail. Then the state has both a 
well defined energy and well defined extent in real space. 
Such a situation certainly exists in a-Si provided that the 
defect is really isolated (in space) from resonant defects. 
From this point of view, the highly localized states near 
the middle of the gap arising from single defects in large 
cells are examples of especially spatially compact chister 
states. As one considers energies closer to the valence and 
conduction bands (in the band tail regions), it becomes 
difficult to compute the cluster states exactly, since the 
density of states increases (which implies that the num- 
ber of cluster states is also greater) so that the probablity 
is large that the tails of resonant defects overlap. If there 
is such an overlap between resonant defects (clusters), 
simple quantum mechanics guarantees that eigenstates 
will consist of mixtures of cluster states with formation 
of "bonding and antibonding" levels. Another way of 
saying it is that we arc in the curious position of: "Given 
the eigenvectors, what is the most localized basis with 
smallest number of elements which can be superposed 
to reproduce the eigenstates?" Such cluster states have 
much in common with our recent work on generalized 
Wannier functions in amorphous materials [10]. 

Our qualitative picture of the local to extended transi- 
tion in a-Si, based on these calculations is the following: 
As severe distortions are rare, clusters stemming from 
such distortions are probably isolated from each other 
and if isolated, are localized energy eigenstates. For less 
severe distortion, the probability of occurrence increases, 
and the size of associated clusters is also larger. Then the 
chance of finding another cluster of similar energy in the 
neighborhood increases. As the distortion becomes less 
severe, eigenstates will consist of mixtures of two, three, 
or more clusters. At some point, clusters can always find 
"overlapping energy partners" and they mix together to 
enable electronic connectivity (extended eigenfunctions). 
This state of affairs can be identified with the "mobility 
edge" . 

With the preceding introduction, it is natural to write 
down a simple "theory" utilizing these ideas. A reason- 
able form for this model is: 

H = Y,Ea\a){a\+ ^ |a)W(/3| (1) 

a a/3(c«#/3) 

where, basis state \a) is a localized cluster {which may 
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involve many atoms). In this representation, the ba- 
sis functions arc localized energy eigenstates with eigen- 
value Ea of the Haniiltonian in the absence of other 
defects with which mixing occurs. is determined 

by the distortion. In real amorphous solids, the clus- 
ter states may significantly overlap, represents the 
coupling between cluster \a) and \(3). In the spirit of 
Hiickel theory, we can take: ^ap ~ {Ea + Efj)Sai3/'2. 
for Safj = {oi\j3). Then, in first order perturbation the- 
ory, the formation of eigenstates of H from these clus- 
ters becomes obvious; the first order correction to the 
zeroth-order (cluster) state \a) is X^/j^^^q ^0/31/3), where 
Tap = {Ea + E/3)S/3a/'2{Ea — Ej3). The strong mixing for 
small energy denominators {Ea « S/3), and the role of 
the overlap is indicated. For reasons that are apparent, 
we have called this the "resonant cluster proliferation" 
model of the Anderson transition [6] . 

Since the cluster basis states la) are not precisely 
known, no calculations have yet been attempted in this 
representation. It is however easy to imagine a calcula- 
tion in which model cluster states are introduced with 
characteristic decay lengths; and the role of the amor- 
phous net is largely encoded into (^afj- Even in the ab- 
sence of new explicit calculations, I assume that the "res- 
onant cluster" model is useful to understanding the An- 
derson transition, and the role of thermal disorder in a 
topologically disordered network, as I describe next. 

4. Thermal modulation of electron states and en- 
ergies 

There is abundant experimental evidence that lattice 

vibrations play an important role in the dynamics of elec- 
trons in amorphous materials. Among other examples, 
Cohen and coworkers [11] observed a pronounced tem- 
perature dependence of the Urbach tails in a-Si:H (the 
conduction tails showing a very strong linear variation 
in decay parameter with temperature). Of course elec- 
trical conductivity is well known to be very temperature 
dependent, and usually has multiple distinct regimes ac- 
cording to different conduction mechanisms [12]. From 
this point of view it is unsurprising that the electron en- 
ergies and states can be very time and temperature de- 
pendent. This point has been independently recognized 
by Arkhipov and Adriaennsens [13] in their studies of 
carrier transport. 

We study thermal modulation of the eigenvalues and 
eigenvectors by tracking their (adiabatic) time develop- 
ment over the course of a few picosecond (lO^^^s) simu- 
lation. As illustrated in Fig. 1 in Ref. [5], there is strong 
time dependence of the LDA eigenvalues in the vicinity 
of the optical gap. The Fermi level is near the middle of 
the gap and several states near the Fermi level are ap- 
propriately described as band tail states. These are much 
like the states which would be responsible for conduction 
in doped a-SiiH. As in earlier work [14] there is a roughly 
linear relation between root mean square (rms) tempo- 
ral fluctuation and temperature. As expected, the higher 



temperature simulation leads to larger excursions in the 
positions of the energy eigenvalues. Note for 300K that 
the Lowest Unoccupied Molecular Orbital (LUMO) fluc- 
tuates in time by about ~0.3eV, much larger than ther- 
mal energies (~10 meV). States deeper into either the 
valence or conduction bands show progressively less ther- 
mal modulation because they are less localized (we have 
noted [15] a very strong correlation between the rms fluc- 
tuation in the energy eigenvalues due to thermal disorder 
and the inverse participation ratio, a simple measure of 
localization in the T=0 model). The localization "ampli- 
fies" the electron-phonon coupling. Also, the conduction 
states fiuctuatc more than the valence states (suggesting 
that the conduction tails arc more sensitive to thermal 
disorder than the valence tails which originate primarily 
from structural disorder), in pleasing agreement with to- 
tal yield photoemission experiments [11] and earlier the- 
ory work [14] . 

Elsewhere [5], we have published color figures depict- 
ing the fluctuations in electron states near the gap for 
the system and dynamics described above. There are 
very substantial changes in the LUMO state in partic- 
ular; There is a clear tendency for the LUMO state to 
alternately "accumulate" on a strained part of the net- 
work, sometimes becoming localized, but also occasion- 
ally developing a substantially more extended "string- 
like" character. These are not the only two recognizable 
structures, but rcc;iu' most frequently. The time between 
"characters" is not predictable, though it is of order tens 
to hundreds of fs. We have posted an animation of this 
state on the world wide web [16]. 

We have argued in the preceding section that struc- 
tural disorder in a-Si gives rise to localized states with 
energies in the band tails. These system eigenstates can 
involve many atoms and can have a Byzantine [17] struc- 
ture. The "simple physics" of this study is that the strong 
(compared to midband electrons) electron-phonon cou- 
pling for localized band tail states is sufficient to cause 
strongly time/temperature dependent quantum mechan- 
ical mixing of cluster states when the thermal disorder is 
"just right" to make their energies nearly degenerate pro- 
vided that they have some overlap in real space. Strong 
mixing of course implies less localization and thus bet- 
ter prospects, at least while the more extended state sur- 
vives, for conductivity and optical transitions. This work 
shows that transport and optical calculations based only 
on T=0 results can be quite misleading. Mott [12] and 
others have made fundamental contributions to the the- 
ory of transport in a noncrystalline medium; for exam- 
ple,variable range hopping. In the kind of simulation we 
present here, I can estimate the conductivity, including 
its temperature and frequency dependence directly from 
the electronic states through an appropriate thermal av- 
erage of the Kubo formula [18]. It is also a complement 
to the phenomenological models of transport [13,19]. In 
the latter work, transport is modeled as a hopping be- 
tween localized tail states. My work can be viewed as 
an explanation of the precise nature of the states among 
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which electrons are hopping (the very comphcated states 
of Ref. 6). The waiting time between hops must be re- 
lated to the time between eigenvalue "close encounters" 
near the Fermi level. It also points at an atomistic level 
to the dynamics of bandtail defects and their kinetics. 

The consequences of this work can be stated another 
way. If \i) (I/)) are initial (final) electronic states with 
energy Ei (Ef), then for an electronic transition in a- 
Si, a Fermi golden rule argument leads quicky to the 
conclusion [12] that the transition rate is proportional 
to \{i\f\f)\'^S{Ef - Ei - fko), where f is a perturba- 
tion inducing the transition (to first approximation a 
momentum operator) and w is the frequency of an ex- 
ternal probe. Both the energies in the S function and 
the transition matrix elements arc affected by the instan- 
taneous details of the structural disorder, and as such 
transition probabilities are also strongly dependent on 
the time and temperature. The consequences of this to 
transport are under investigation; the discussion here is 
based upon first-order time dependent perturbation the- 
ory, which for the very strong electron-phonon coupling 
we discuss, could be inadequate. 

5. Heuristics 

Earlier work has shown that it is useul to link the ther- 
mal fluctuation of the LDA energy eigenvalues near the 
band tails to the extent (localization) of the band tails 
in amorphous Si (as separately measured in total yield 
photoemission experiments). Here I give a simple model 
which makes this connection between fluctuating electron 
energies and lattice vibrations a bit more explicit. 

Consider a particular electronic eigenvalue, A„, say in 
one of the band tails in a-Si. To estimate the sensitiv- 
ity of A„ to a coordinate distortion (supposedly ther- 
mally induced), we can use the Hellmann-Feynman the- 
orem [20], which gives 9A„/9Rq = {ijjn\dtl / dlla\ipn) 
(for this to be valid, I must assume that the basis is 
flxed (not moving with the atoms) and that the \tpn) 
are exact eigenvectors of H; see Ref. [8] for a more gen- 
eral case). Then clearly for small distortions R^, we 
have 6Xn ^ ^^(^/'„|9H/9Ra|^„)(5Ra. Here, R is the 
3N vector of displacements for all of the atomic coor- 
dinates from equilibrium. If the displacements SRait) 
arise from classical vibrations, then one can also write 
5Ra{t) = J2uj MT,u}) cos[Lut + (pujxaico), whcrc 
CO indexes the normal mode frequencies, A{T,w) is the 
temperature dependent amplitude of the mode with fre- 
quency u), (poj is an arbitrary phase, and Xait-o) is a 
normal mode with frequency iv and vibrational displac- 
ment index a. Using a temperature dependent ampli- 
tude A^{T,uj) = kBT/2Muj'^, it is easy to sec that 
the trajectory (long time) average of the expression for 
SXl is: {SXD ~ (fcB^/4A^)EJ5(a')/a;)^ where 
= E„(V-'«|5H/aRa|'(/;n)Xa(t^)]. It would be 
straightforward for a particular collection of vibrational 
states to include the correct Bosc terms to obtain a result 
valid at low temperatures T < Qd (for a salient De- 



bye temperature). These expressions give a transparent 
expression for the thermally induced electron modula- 
tion as driven by the lattice-electron coupling S(a;). If 
we follow Ref. [14] and make the coarse approximation 
of equating ((5A^)^/^ to the Urbach decay parameter [11], 
then this model would predict a square root dependence 
of the decay parameter with temperature, perhaps not 
easily distinguished from the linear dependence reported 
in Ref. [11]. 

It is easy to see why there is a correlation between 
an eigenstate's localization (as measured for example by 
inverse participation ratio) and its thermal rms fluctu- 
ation [15]. Note that S(w) will in general consist of a 
sum of many terms (for different a) ; the individual terms 
have no preferred sign, so that adding a large number 
of terms of comparable magnitude will lead to cancel- 
lation and a small sum. On the other hand, if only a 
very small number of terms are nonzero (as for the case 
when the electron state |f/'n) is well localized), then there 
will be less cancellation and a larger contribution to the 
sum (and therefore to the fluctuation of the eigenvalue). 
This model is limited in many ways (it is classical, as all 
the modes are equally populated from classical cquipar- 
tition), it is obviously strictly harmonic and we arc as- 
suming no electronic level crossings or other departures 
from adiabatic (Born-Oppenheimer) dynamics [5]), but 
captures some of the right temperature dependent elec- 
tronic effects observed. 

The thermal effects on the electron states are more dif- 
ficult to calculate with this approach, mostly because of 
problems with degeneracy if one tries to use perturbation 
theory. Still, the underlying "physics" is quite simple, 
just the effects of resonant mixing for close approach (in 
energy) of eigenvalues originating in cluster states with 
some overlap in space. In the language of the model 
Hamiltonian above, this means that there can be large 
mixing of cluster basis states when the energy denomi- 
nator becomes small; here the thermal fluctations can in- 
duce small energy denominators (and therefore mixing). 
Of course a mixed state involves more cluster states than 
a pure cluster state, so that it is more spatially extended. 

5. Post adiabatic atomic dynamics 

The calculations and discussion of the preceding sec- 
tions has been based on the adiabatic approximation 
(that the electrons are always in their instantaneous 
ground state determined by the external potential (eg 
the atomic nuclei), and that the forces "felt" by the nu- 
clei arose from the instantaneous electronic structure of 
the system. 

To see that this might be an imperfect assumption, 
consider a doped system with a small concentration of 
dopants for which the Fermi level is pushed into one 
one of the band tails. Then the thermal modulation 
of the atomic coordinates may produce level crossings 
or close approaches of the LDA energies. In the doped 
case, there would normally be (occupied) levels below 
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the Fermi energy and (unoccupied) levels above. In the 
case of a level crossing at the Fermi level, the adiabatic 
approximation would immediately transfer charge from 
one state to another (even if spatially remote), which 
would induce changes in the dynamics (since the lev- 
els with changed occupation would typically be at least 
somewhat localized); and perhaps could lead to a struc- 
tural change, even a long lived change if the new struc- 
ture was "self-trapping" . In practical calculations such 
charge transfer would not have be quite so large an affect, 
since one would usually smear the Fermi function slightly 
from the T = step form. The dynamics are then af- 
fected by such level crossings or close approaches even 
in the adiabatic approximation. A better theory [21] is 
due to Allen, in which interatomic forces are computed 
from mixed state wave functions obtained by a direct 
integration of the time- dependent Schrodinger equation 
(rather than the "pure" states computed anew at each 
new time step for the atomic dynamics in the adiabatic 
picture) . This approach was formulated to model the dy- 
namical response of a system to light (in which the sim- 
plest model of light-solid interactions would just involve 
promotions of electrons to low- lying conduction states), 
but the formalism is similar (and much simpler) for the 
case of purely phonon-induced level crossings as in our 
problem. It is possible in fact that the adiabatic ap- 
proximation may often be satisfactory for this kind of 
problem, but the question has not yet been properly in- 
vestigated because of difficulties with post-adiabatic cal- 
culations. One case for which nonadiabatic dynamics is 
essential is in the modeling of the thermalization of ex- 
cited carriers to tail states. This process has never been 
realistically modeled despite its manifest importance to a 
host of problems in crystalline and non-crystalline semi- 
conductor physics. More sophisticated methods for non- 
adiabatic dynamics are discussed in the quantum chem- 
istry literature [22] ; while fundamentally sound, these are 
currently too difficult to implement for the large model 
systems we must consider here. 

6. Conclusions 

We have used current computational techniques to 
demonstrate the nature of the localized to extended tran- 
sition in amorphous Si, and then took the additional 
step of adding thermal disorder and found that local- 
ized states could be time dependent. A simple analytic 
model was then presented, which explains some of the 
general temperature-dependent features observed in ex- 
periments. Finally we speculated about the importance 
of non-adiabatic dynamics to these calculations. 

I note in closing that the effects reported here are pre- 
sumably not peculiar to a-Si; they can be anticipated for 
any disordered insulator. 
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